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Global maps (l.e.» I.R. emitted flux, net flux, and albedo) of 
processed NIMBUS-6 data for July, 1977. Units for fluxes 
(watts/m^), for albedo (percent). 

Same as 1, except for August, 1977. 

Same as 1, except for September, 1977. 

Same as 1, except for October, 1977. 

Same as 1, except for November, 1977. 

Same as 1, except for December, 1977. 

Same as 1, except for January, 1978. 

Same as 1, except for February, 1978. 

Same as 1, except for March, 1978. 

Same as 1, except for April, 1978. 

Same as 1, except for May, 1978. 

Scale and geography for global maps (I.e., Figs 1-11) 

Zonal averages of net radiation flux for April. 

Same as 1-12, except for May. 

Same as 1-12, except for June. 
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Zonal mean emitted flux (watts/m ) for Dec - Jan period 
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Zonal mean net flux (watts/m ) for Dec - Jan period 

Net solar radiation at the surface and top of the model’s atmosphere 
for January 

Observed and calculated mean zonal wind 500 mb 
Observed (200 mb) and calculated (250 mb) mean zonal wind 
Observed (700 mb) and calculated (750 mb) meridional wind 
Observed and calculated pressure velocity u 500 mb 
Calculated precipitation and mixing ratio day 91-120 




During the first seven months of research on Radiation Budget Measure- 
ment/Modeling Interface we have laid plans and completed early work In three 
areas: (a) Improved Earth Radiation Budget Data Sets, (b) Numerical Model 

Experiment Definition and (c) Review of Concurrent Research by Other Scientists. 
All of these areas of work will be reported In detail via special reports now In 
preparation. Figures Included In the following sections of the present report 
Illustrate the work In prog 'ess. 


1.0 ImDroi^ad Earth Radiation Budqat Data Sat« 


As a guide to selection of climate modeling experiments we are continuing 
our early work by processing and analyzing the NIMBUS-6 data. The first two 
years or NIMBUS-6 data as well as earlier data« which have been previously 
processed, were reported by Stephens ^ ai. (1981). In recent months we com- 
pleted the processing of the third year of NIMBUS-6 data. Unfortunately, the 
data quality from June, 1978 ough October, 1978 was Insufficient to provide a 
meaningful analysis of the earth's radiation budget for this period. Thus, the 
third year of data extends from July, 1977 to May, 1978. The processed results 
for the "third year" of the NIMBUS-6 data are shown In Figures 1-1 through 1-11. 
The scale and geography for these global maps Is shown In Figure 12. These data 
are being combined with earlier data (see Table 1) to form the most up-to-date 
climatology of the earth's radiation budget from which the model work noted 
below will be guided. Special attention has been directed to the April - June 
period, which Is hypothesized by the P. 1. to be a principal time of potential 
Interannual variability In the radiation budget. 

From the processed data generated from the satellites referenced In Table 
1, we have computed the zonal averages of net radiation for the three months 
mentioned above (Figures 1-13 through 1-15). It Is still uncertain whether the 
differences observed In these zonal averages are. In fact, manifestations of an 
Interannual variability within the earth's radiation budget, or If they are 
merely artifacts of the various satellite sampling systems. 

Assuming satellite problems to be minimal, the Northern Hemisphere winters 
of 1975/76, 1976/77 and now 1977/78 are particularly Interesting. In contrast 
to Campbell (1981) who surmised only minor differences between 75/76 and 76/77, 
our preliminary results from this research show 75/76 and 77/78 to be strikingly 
similar In terms of radiation budget - with 76/77, the winter of major North 
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Figure 1-4 
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Figure 1-6 
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Figure 1-7 
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Figure 1-8 
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Table 1. Cliroitologlcal list of earth orbiting satellites from which present radiation measurements 
were taken. The approximate local time at whlcli cacli satellite crossed the equator during 
daylight hours appears In parenthesis. EX - ' xperlniental, N2 - Nimbus 2, N3 -> Nimbus 3, 

N6 - Nimbus 6, E 3 - E.sf-'J and E 7 • fssa 7. 
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LRTITUDE (DEGREES) 




American weather anomalies, distinctly different from the preceeding and follow- 
ing years. Differences can be noted In the December-January zonal means for 
emitted and net fluxes for these years (Figures 1-16 and 17, respectively). The 
fact that the emitted flux (Figure 1-16) In Dec. 76 - Jan. 77 was larger at 
higher latitudes (60°N - 90°N) and lower at mid-latitudes (30®N - 90®)N In 
comparison to the other years may be Indicative of the strong and persistent 
atmospheric "blocking" patterns which character Ized this period In the Nothern 
Hemisphere. For under such blocking conditions, the Increased meridional cir- 
culation would allow anomalous amounts of warm air to be advected Into higher 
latitudes which. In turn, would result in Increased longwave emission. Conver- 
sely, at mid-latitudes, the Increase In cold air advectlon associated with the 
blocking flow would result In decreased longwave emission. Since the NCAR 
Community Climate Model (XM) has recently replicated the blocking situation 
found over North America during the 1976/77 winter (Blackman, 1981), our model 
experiments planned for the remaining portion of the research period of this 
contract have additional Impetus potential. 

A more complete documentation of the new NIMBUS-6 results, as well as 
preliminary discussion. Is In preparation (CleslelskI, Campbell, and Yonder 
Hear, 1981). This special report under the present contract will be patterned 
after Ellis and Yonder Haar (19/6) and Campbell and Yonder Hear (1980). 

2.0 Numerical Model Exp eriment Definition and Tests 

Two complementary models are presently under development. The first Is 
a statistical -dynamical model patterned after the recent work of Ashe (1979). 

It contains very crude vertical resolution (only two levels; but In principle, 
more can be added) and an arbitrary horizontal resolution. The horizontal 
dependence Is represented by spectral components, which have recently been 
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demonstrated superior to a finite difference grid. Because the dissipation Is 
often modeled by horizontal diffusion. It Is useful to use spherical harmonics 
ar- the spectral basis functions. 

In this first model approach, the diabatic forcing Is assumed temporally 
constant (e.g., a seasonal average) and we solve for the steady circulation as a 
response. Transient motions (I.e., weather) are parameterized as dissipative to 
tho steady circulation. To make the model tractable, the results understand- 
able, and the energetics simple, a dynamical approximation to the primitive 
equations Is assumed. At present the "linear balance" system Is being used 
because It Is the simplest quasl-geostrophic-type system which Is valid 
globally. Despite the "I'near balance" label, the vorticity and thermodynamic 
budgets are fully non-l Inear. This Is the chief distinction from the generation 
of linear models (e.g., Sankar-Rao and Saltzmann) used In the 1960's; In those 
earlier studies, the mean zonal flow was specified a priori Independent of the 
standing eddies. Hence, a large part of the general circulation was assumed , 

and not necessarily consistent relative to the eddies. Here we take the much 

« 

more satisfying approach of finding both the mean flow and the standing eddies 
which result from a prescribed stationary forcing, consisting primarily of 
radiative heating/cool I ng, latent h'^at release, a.''d sensible/ latent heat flux 
from the surface. 

Unfortunately, this Important aspect of reality also makes the problem much 
more difficult to solve, because a non-l Inear system of equations mucT be 
solved. There Is no universally applicable method for finding the non- 1 1 near 
solution; Indeed, there Is not even a guarantee of a unique solution or of any 
solution at all. Because of Ashe's results, we are confident that a solution 
can be obtained. At present, we are attempting three methods of solution: 
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(1) Iterate In time from a simple Initial condition to the asymptotic 
steady state. 

(2) Use a standard IMSL non-linear equation solver. (Unfortunate I y» this 
method depends crucially on Initial conditions, »hlch oay not be 
sufficiently precise.) 

(3) First calculate a linear response and then use this solution as the 
Initial guess In the Iterative methods (1) and/or (2). 

We still have a significant amount of work to do In this area In order to 
obtain the non-l Inear solution. Because of the notorious difficulties with 
non-l Inear systems of equations, there Is no guarantee of success. 

The second model approach makes similar dynamic approximations, but 
explicitly calculates the time evolution of the flow. Thus, time-dependent 
radiation can be taken from observations (Including the seasonal cycle and the 
Interannual variation), applied as a forcing, and the response of the atmosphere 
computed. The Individual "synoptic" flows can then be averaged to obtain 
monthly means, seasonal means, and Interannual variations. Thus, we will be 
able to find whether or not the atmospheric circulation Is sensitive to exter- 
nally specified changes In radiative forcing. In turn, this will guide our 
understanding and Interpretation of the present and future radiation budget 
measurements from satellites. 

The prel li.'«lnary developrent of the latter model Is being done by Adel 
Hanna, a Ph. D. student, utiuer the auspices of other funding sources. Duane 
Stevens Is advising him In this research effort. As this model Is In an early 
development stage, results are of a very preliminary nature. Further "tuning" 
of the physical parameterization will probably be required In order to attain an 
approximately veritable simulation of the annual average and seasonal I y-vary I ng 
circulation. The following represents some of these preliminary results. 
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Th« second modol under discuss I or Is a two- 1 eve I global "linear bo lance" 
node I • In the horizontal domain a rhcmboldal spectral truncation Is assumed 
(presently, we truncate at zonal wave no. 9). In this section the results of 
time Integration of the model up to 120 days, to simulate the Ja-^uary circula- 
tion, are discussed. 

The model was Initiated assuming an atmosphere at rest and with constant 
moist adiabatic lapse rate (6.5^/km). Assuming a perpetual January forcing, the 
solar Inclination Is fixed to that of the first of January. Figure 2-1 shows 
the net solar radiation at the top of the model's atmosphere and the net solar 
radiation absorbed by the earth's 'surface at different latitudes In the model. 
Except the Artic Ocean, sea surface temperature were prescribed as the January 
climatological values. Orographic of feds are parameterized through the verti- 
cal motion at the 1000-mb level. 

The winter observational estimations for different variables were taken 
from Newell jbI aI* (1971). Data for the average January were taken from Oort 
and Rasmusson (1971). 

a) The Average Zonal Mind 

The zonal wind component at different grid points was averaged In 
latitude and tltre (30 days). In 'fhe northern hemisphere (winter hemisphere) 
the average 300-mb zonal wind compared well with tha observed pattern (Figure 
2-2). The same for the 250-mb distribution (Figure 2-3). The observations 
show a westerly Jet at 30^ (20 m sec~^ at 300 mb and 38 m sec**^ at 230 mb). 

The calculated values place the maximum zonal wind at 23^ with maxima of 13 m 

sec~^ and 23 m sec~^ for the 300-mb and 230-mb cases, respectively. It can also 
be noticed that the calculated averages for the period 61-90 days show nearly 
the SMie values as the averages for the periods 91-120 days In the case of 300 

mb, also the same agreement between the 31-60 days average and 91-120 days In 
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the case of the 250 mb. This Indicates that the zonal momentum Is reaching a 
statistical steady state. 

In the southern hemisphere (summer hemisphere) the calculated values have 
less agreement with observations than the northern hemisphere case. In both the 
calculated 500-mb and 250-mb cases, the westerly Jet In the extratropical 
southern hemisphere Is much damped and shifted. On the contrary, the equatorial 
easterlies aro amplified and broadened to reach the southern hemisphere mid- 
latitudes. The variability between different time averages for both 500 mb and 
250 mb may suggest that values at this region need a longer simulation than 
already have been done. 

b) Meridional Wind Component 

The zonal I y averaged meridional wind component Is being calculated from 
the zonal ly-averaged velocity potential. The latter Is at least one order of 
magnitude smaller than the velocity stream function (In middle and high 
latitudes). Even observational estimates show a large variability between 
winter and January cases (Figure 2-4). In the northern hemisphere the calcu- 
lated averaged meridional velocity at 750 mb changes Its sign corresponding to 
the vertical mass flux at 500 mb (Figure 2-4), validating the mass continuity. 

In the southern hemisphere the areas of equatorward flux (positive v) are not 
simulated since the vertical mass flux Is mainly upward (negative w). 

c) Vertical Velocity 

-4 

The observed and calculated 500 mb vertical velocity field (units 10 

-1 -5 -1 

mb sec , 10 mb sec , respectively) Is shown In Figure 2-5. Again In the 
northern hemisphere the calculated values resemble the observed phase but with 
smaller amplitudes. The midlatitude oscending motions and the subtropical 
subsidence are well fitted with the observations. The calculated subtropical 
subsidence occur*; In a rather broad latitude band In comparison to that which Is 
observed. This featire may account for the same mass flux to counter the 
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Fig. 2-4 0bserved(700 mb) and calculated (750 mean meridional wind 





decrease of amplitude. In the southern hemisphere maximum shifting with the 
ITCZ Is shifted to 30^S with a very weak subsidence In the extratropical 
latitudes. 

d) Moisture Budget 

The model uses a simple moisture budget equation In the lower layer 
with the assumption that the upper -ayer Is dry. Precipitation occurs with 
release of latent heat of condensation when the relative humidity of the lower 
layer exceeds 80$. The model's atmosphere Is convective I y adjusted If the 
temperature lapse rate, after the release of latent heat, exceeds the saturated 
adiabatic lapse rate. Figure 2-6 shows the zonal ly averaged precipitation and 
mixing ratio (gm/gm); both have a reasonable distribution except the subtropical 
maxima at 30^. This Is a result of the critical relative humidity assumed for 
precipitation (80j^). At the subtropical latitudes the relative humidity may 
exceed 80|^ without any precipitation. 



3.0 Rflviaw of Coneurrant Radlwtlon Budqat Raaaarch by Qthar Sclantlst^ 


During the last year new papers and reports by Ohring, 6. Smith fit 
Gruber, Hansen fit fit., and others have continued the study of earth radiation 
budget and related modeling. We are reviewing these papers as they relate to 
our work on this project. They will be referenced In our reports In prepara- 
tion. 
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